RANDOM WALKS ESTIMATE LAND VALUE IN 
CITIES 

Ph. Blanchard &c D. Volchenkov 

Bielefeld - Bonn Stochastic Center, University of Bielefeld, Postfach 100131, D-33501, 
Bielefeld, Germany 

E-mail: blanchardSphysik . uni-bielef eld . de , 
volchenkOphysik . uni-bielef eld . de 

Abstract. 

Expected urban population doubling calls for a compelling theory of the city. 
Random walks and diffusions defined on spatial city graphs spot hidden areas of 
geographical isolation in the urban landscape going downhill. First-passage time to a 
place correlates with assessed value of land in that. The method accounting the average 
number of random turns at junctions on the way to reach any particular place in the 
city from various starting points could be used to identify isolated neighborhoods in 
big cities with a complex web of roads, walkways and public transport systems. 
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1. The problem of massive urbanization 

Multiple increases in urban population that had occurred in Europe at the 
beginning of the 20 th century made urban agglomerations suffer from the problems 
of urban decay such as wide-spread poverty, high unemployment, and rapid changes 
in the racial composition of neighborhoods. Riots and social revolutions happened in 
response to conditions of urban decay in many European countries established regimes 
affecting immigrants and certain population groups de facto alleviating the burden of 
the haphazard urbanization by increasing its deadly price. 

More than half of world's population, 3.3 billion people, is now living in cities, and 
the figure is about to double by 2030 [Tj. Urban sprawl in US covered 41, 000 km 2 in 20 
years, the area that equals that of the state of Switzerland [2]. In Europe, urban sprawl 
has covered 8, 000 km 2 in 10 years the area that equals the territory of Luxembourg [3]. 
Unsustainable pressure on resources causes the increasing loss of fertile lands through 
degradation and climate changes through the increasing carbon emissions warming the 
earth's atmosphere. City development planners will face great challenges in preventing 
cities from unlimited expansion. 
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Global poverty is in flight becoming a primarily urban phenomenon in the 
developing world. 2 billion new urban settlers in the next 20 years will live in slums on no 
more than $1 a day, adding to 1 billion already living there [3]. Faults in urban planning, 
poverty, redlining, immigration restrictions and clustering of minorities dispersed over 
the spatially isolated pockets of streets trigger urban decay, a process by which a city 
falls into a state of disrepair. The speed and scale of urban growth in that require urgent 
global actions to help cities prepare for growth and to avoid them of being the future 
epicenters of poverty and human suffering. 

2. The scope of study and conclusions 

In the present paper, we investigate spatial graphs of several compact urban patterns: 
the city canal networks of Venice and Amsterdam, as well as the almost regular array 
of streets in Manhattan. 

In the next section (Sec. [3]), we briefly discuss the principles which is are used 
for encoding urban environments into geometrical graphs. In Sec. HI we develop the 
theory of linear automorphisms for undirected graphs and demonstrate that they can 
be interpreted as random walks on graphs. By the way, many geometrical properties 
of graph may acquire probabilistic interpretations. In particular, we are interested in 
first-passage times to the nodes of graphs - the expected number of steps required for 
random walkers starting from an arbitrary node of the graph to reach the given node 
for the first time. The first-passage time is the important characteristic of the node 
quantifying its isolation in the complex graph structure. 

Sociologists think that isolation worsens an area's economic prospects by reducing 
opportunities for commerce, and engenders a sense of isolation in inhabitants, both 
of which can fuel poverty and crime. It is well known that many social variables 
demonstrate striking spatial distribution patterns, and therefore may be detected and 
predicted by a structural analysis. In the last Sec. El we discuss two examples of such a 
type: the spatial isolation of Venetian Ghetto, the small area of Venice in which Jewish 
people were compelled to live under the Venetian Republic, and the value of land in 
the urban pattern of Manhattan. In particularly, we compare the tax assessment rates 
of places in Manhattan with first-passage times to them. The data positively relate 
the geographic accessibility of places in Manhattan with their 'unearned increments' 
quantifying the level of appreciation in value of the site. 

The proposed method could be used to identify isolated neighborhoods in big cities 
with a complex web of roads, walkways and public transport systems. 

3. Spatial networks of urban environments 

In traditional urban researches, the dynamics of an urban pattern come from the 
landmasses, the physical aggregates of buildings delivering place for people and their 
activity. The relationships between certain components of the urban texture are often 
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measured along streets and routes considered as edges of a planar graph, while the 
traffic end points and street junctions are treated as nodes. Such a primary graph 
representation of urban networks is grounded on relations between junctions through 
the segments of streets. The usual city map based on Euclidean geometry can be 
considered as an example of primary city graphs. 

In space syntax theory (see [HE]), built environments are treated as systems of 
spaces of vision subjected to a configuration analysis. Being irrelevant to the physical 
distances, spatial graphs representing the urban environments are removed from the 
physical space. It has been demonstrated in multiple experiments that spatial perception 
shapes peoples understanding of how a place is organized and eventually determines the 
pattern of local movement, [6]. The aim of the space syntax study is to estimate the 
relative proximity between different locations and to associate these distances to the 
densities of human activity along the links connecting them, EJ |9] . The surprising 
accuracy of predictions of human behavior in cities based on the purely topological 
analysis of different urban street layouts within the space syntax approach attracts 
meticulous attention [TO] . 

The decomposition of urban spatial networks into the complete sets of intersecting 
open spaces can be based on a number of different principles. In [11], while identifying 
a street over a plurality of routes on a city map, the named-street approach has been 
used, in which two different arcs of the primary city network were assigned to the same 
identification number (ID) provided they share the same street name. 

In our paper, we take a "named-streets" -oriented point of view on the decomposition 
of urban spatial networks into the complete sets of intersecting open spaces following 
our previous works [321 H3] . Being interested in the statistics of random walks defined 
on spatial networks of urban patterns, we assign an individual street ID code to each 
continuous segment of a street. The spatial graph of urban environment is then 
constructed by mapping all edges (segments of streets) of the city map shared the 
same street ID into nodes and all intersections among each pair of edges of the primary 
graph into the edges of the secondary graph connecting the corresponding nodes. 

4. Theory of linear automorphisms for undirected graphs 

4-1. Graphs and their linear automorphisms 

Consider a transitive permutation group V(V) defined on a finite ordered set of 
objects V, \V\ = N. Its representation consists of all N x N orthogonal matrices Yl a 
such that (n o .) i ai = 1, and (II CT ) i . = if j ^ oi. We can define an induced action of 



Given a binary relation of adjacency "~" defined on the finite set V, we denote its graph 
by G and conveniently represent by its adjacency matrix A - the N x N matrix where 




v, w e V. 



(1) 
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the non-diagonal entry Aj = 1 if and only if % ~ j, and is zero otherwise. 

Among all possible permutation groups V(V) defined over the set V, there is one 
'compatible' with the binary relation "~" that is the automorphism group of the graph 
Aut(G) such that for any o G Aut(G), uo ~ vo~, if and only if u ~ v. A permutation 
a G Sym(G) belongs to Aut(G) if and only if its matrix representation 11^ commutes 
with the adjacency matrix of G [T4"] . 

[A, n CT ] = 0. (2) 

A linear function of the graph G defined on its adjacency matrix, 

JV 

/(A);. = F ijtSl A sl , F ihsl g R, (3) 

s,l=l 

belongs to the representation of the automorphism group Aut(G) in the group of linear 
transformations ([3]) if 

nj / (A) n a = f (nj a n a ) & £ 4* = E / ( A h » ( 4 ) 

for any a G Aut(G). In [15] . we have shown that any function defined on a simple 
undirected graph G and satisfying (j3J) can be associated to a stochastic process, 

/ (A) y = (1 - /3) h + —Aj, «» = J] ^i, (5) 

where /3 G [0, 1] is an arbitrary parameter, and Ki is the degree of the node i G V. The 
operator / (A) is the transition operator of "lazy" random walks for which a random 
walker stays in the initial vertex with probability 1 — 0, while it moves to another node 
randomly chosen among the nearest neighbors with probability /3/«j. In particular, for 
(3 = 1, the operator / (A) describes the standard random walks extensively studied in 
classical surveys [16]. [TT]. 

4-2. Algebraic geometry of linear automorphisms 

The attractiveness of random walks methods relies on the fact that the distribution 
of the current node of any undirected non-bipartite graph after t > 1 steps tends to a 
well-defined stationary distribution, 

TTi = ^ ~ — , (6) 

which is uniform if the graph is regular. The stationary distribution of random walks 
OH]) defines a unique measure on the set of nodes V with respect to which the transition 
operator ((J3J) for /3 = 1) is self-adjoint, 

T = 1 -(n^fn-^ + n-^fn^), (7) 

where f T is the adjoint operator, and n is defined as the diagonal matrix n = 
diag(7ri,...,7Tjv)- 
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In the canonical basis associated to the nodes, {(0, . . . , lj, . . . , 0)} i=1 , with 1 at the 
i-th position, for every node i 6 G, the self-adjoint operator ([7]) is represented by a real 
positive stochastic matrix. It follows from the Perron- Frobenius theorem [18] that its 
maximal eigenvalue is simple and equals fi\ = 1, and all components of the eigenvector 
ipi belonging to the largest eigenvalue fix are positive, ipi,i — y/^U > 0, Vi G V. All 
eigenvalues are real 1 = H\ > \ii > ... > /xjv > — 1, with orthonormal eigenvectors 

N N 

E <i = E << = L ( 8 ) 

1=1 fc=l 

forming an orthonormal basis in Hilbert space 1-L(y). For eigenvalues of algebraic 
multiplicity r > 1, a number of linearly independent orthonormal eigenvectors can 
be chosen to span the associated eigenspace. 

Eigenvectors map the nodes of the graph onto the (N — l)-dimensional unit hyper- 
sphere, ipf. : V — > S[ N . Let us denote the matrix in which the eigenvectors ipk are 
columns by It is clear that ^ G O(N) where O(N) is the group of all orthogonal 
N x iV-matrices since ^ ,T ^' = 1 and det ^ = ±1. The orthogonal matrix has 
the skew-symmetric property, ^ n ,m = —^m,n, and therefore contains just N(N — l)/2 
independent entries which are determined by iV — 1 autonomous angular coordinates Q\ , 
I = 1, . . . N — 1. If we denote the degree of the node % G V by «j = 2\E\ipf i( then the 
angular coordinates can be defined by 

1 1 

tan 2 6 l+1 = — — E K ^ l=l,---N-2. (9) 

Then the entries of the first eigenvector are given by the standard relations between the 
Cartesian and the hyper-spherical coordinates on the unit hyper-sphere, 

V>i,i = cos 0i, 

/ a -n N ~ 2 ■ a (^) 

V>i,jv-i = cos0iv_i H s=1 smt) S) 

i>i,N = Ils=i siri Qs- 

Since all components of the first eigenvector ip\ = ^/tt are positive, we can construct the 

stereographic projection from the unit hyper-sphere S^ -1 onto the hyper-plane PM.^ , 

with the use of the point ipi = (^i.i, • • • ^i,n) as the center of projection. The projective 

plane PM.^ ^ is then spanned by the basis vectors, 

ft = {^il, k = 2,...,N, i eV, (11) 



so that any vector 7r _1 / 2 v G PM.^ can be expanded into 7r _1 / 2 v = ^2^ =2 VPk where 

^ = l^)(^l (12) 

is the projection operator onto the space associated to the basis vector ip' k . The set of all 
isolated vertices *P of the graph G(V, E) for which 7Tsp = play the role of the plane at 
infinity, away from which we can use the basis defined by fllip as the Cartesian system. 
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4-3. Linear automorphisms and probabilities 

The first eigenvector belonging to the largest eigenvalue //i = 1 of (JZJ) is nothing 
else as the Perron-Frobenius eigenvector which determines the stationary distribution 
of random walks 1161. 



l/>lT = Ipi, = TTi, (13) 

and the Euclidean norm in the orthogonal complement of ipi, Yl^^^Psi = 1 — tTi, is 
nothing else as the probability that a random walker is not found in i. In a continuous 
time jump process, the transition time r is a discrete random variable, with mean 1, 
distributed with respect to the Poisson distribution Po(t). The transition probability 
for the continuous time jump process is 



(14) 



Therefore, eigenvectors from the orthogonal complement of ip\ determine a relaxation 
process toward the stationary distribution of random walks tt described by the set of 
characteristic time scale of the transient process, 

r fe = (l-M fc r\ k>2. (15) 

For each n — 0, . . . , N, we can construct a new vector space /\ n S^ -1 , where the symbol 
A denotes the standard wedge product of vectors in S^ 1 . It is clear that /\° S^ 1 = R, 
A/ Sf -1 = Sf-\ and A fe S? ~ l (2<k<N) consists of all sums £ a a x A a 2 A . . . A a k 
where a e R, a { € 

Given an ordered set of indices X = k 2 , . . . , k n }, 1 < k\ < k 2 < ■ ■ ■ < k n < N , 
the forms 

= ^A.-.A ip kn (16) 
over all possible [ | ordered sets X is an orthonormal basis in /^R^. As usual, 

V n J 

the determinant (|T6|) can be considered as a linear and signed version of hyper-volume: 
the absolute value of the determinant of real vectors is equal to the volume of the 
parallelepiped spanned by those vectors. It is clear that |det = 1, and moreover the 
natural normalization condition holds, 

Z>2 = 1, (17) 



in which the sum is taken over all possible ordered sets of indices X. The 




squared hyper-volumes ip^ can therefore be interpreted as probabilities to find a 
random walker contributing to the relaxation processes determined by the eigenmodes 
X = {h, k 2 ,..., k n }, 

Pr[k u k 2 ,...,k n ] = Vf. (18) 
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The obvious and simplest example is the probability to observe a random walker in 
v during infinite time, Pr[l]| y = ty\ v — ir v , ^2 v tt v = 1> given by the stationary 
distribution of random walks over nodes. 

4-4- Linear automorphisms and first hitting time 

In graph theory, the distance between two vertices in a connected graph is the 
number of edges in a shortest path connecting them. It is also known as the geodesic 
distance because it is the length of the graph geodesic between those two vertices. 
However, given a random walk defined on the undirected connected graph G(V,E), 
each path in that can be characterized by a certain probability to be followed by random 
walkers, so that the shortest path between two nodes may be different from the most 
probable one. In fact, while travelling from the vertex i G V to the vertex j G V, random 
walkers follow with different probabilities all possible paths that connect them. Hence, 
while discussing linear automorphisms of undirected connected graphs, it seems natural 
to generalize the notion of distance between two vertices in the graph in a probabilistic 
sense, namely as the expected number of random steps required to random walkers 
starting at i G V in order to reach j G V for the first time, 

D(i,j) = Uj. (19) 

The above defined quantity is called the first hitting time, [T6] . 

The standard computation of the first hitting time which can be found in [16] is 
based on the fact that while travelling from i G V to j G V, the first step takes us to a 
neighbor v G V of i and then we have to reach j from there, 

Uj = 1 + V] t vj , (20) 

together with the obvious requirement tu = 0. 

Alternatively, we note that random walks defined on finite undirected graphs can be 
naturally related to a diffusion process which describes the dynamics of a large number 
of random walkers. The symmetric diffusion process corresponding to the self-adjoint 
transition operator of random walks ([2J) determines the time evolution of the normalized 
expected number of random walkers, TT- 1 / 2 !!^) G V X N, 

h = Ln, L = 1 -T (21) 

where L is the normalized Laplace operator defined on G. Eigenvalues of L and T are 
simply related by Xk = 1 — fik, k — 1, . . . , N, while eigenvectors of both operators are 
identical. The analysis of spectral properties of the operator ( 121]) is widely used in the 
spectral graph theory, [121 [20] • 

The expected first hitting time can formally be calculated as 
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if the inverse Laplace operator l/L = ^ n>l T n exists. The r.h.s. of ( )22l) is nothing 
else but a convolution of the projection operator onto the difference direction vector P^- 
with the Fredholm kernel of the Laplace equation. 




Figure 1. Schematic representation of a projective action of the inverse Laplace 
operator, : Sf" 1 - {^i} ->• PM^" 1 . 



Although the Laplace operator L has a trivial eigenvalue Ai = 1 — /xi = and 
therefore is not invertible over S^" 1 , it is invertible over S^ 1 — {ipi}, the orthogonal 
complement of the first eigenvector tpi (belonging to the largest eigenvalue of ([7]) 
\L\ = 1) that corresponds to the transient process of random walks toward the stationary 
distribution tt. The orthogonal complement is homeomorphic to the projective hyper- 
plane PM^ -1 ^ constructed by linearly mapping points of the unit hyper-sphere S^' 1 
from ipi as the center of projection. The inverse Laplace operator defined by the kernel 

N 



u 



f— f 1 - Uk 



k=2 ^ k 



in which Vk is the projector (JT2]h is itself a projection operator L 1 : 1 — {^1} — > 
FM W_1 where the k th coordinate of the projective manifold PM^ -1 is subjected to 
the dilatation r k = (1 - // fc )~\ k = 2,...,N (see Fig. Kty . The kernel (E3J) defines 
the Green function (or the Fredholm kernel) describing long-range interactions between 
eigenmodes of the diffusion process induced by the graph structure. 

The convolution with Green's function gives solutions to inhomogeneous Laplace 
equations. It is remarkable that such a convolution with the projection operator 
gives for the first hitting time the expression 

ta = £ tjlz^ktk , (24) 

k=2 1 ^ 

which coincides with the result of |16|. 
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It is well known [T7J US] that the matrix of first hitting times is not symmetric, 
even for a regular graph, and therefore in general D(i,j) ^ D(j,i). However, one can 
consider its symmetrized analog, 

N M 2 _Jj> 2 

d(ij) = D(i,j) + D(j,i) = — , (25) 

known as commute time |16j . the expected number of steps required for a random walker 
starting at i to visit j for the first time and then to return back to i for the first time. 
It is easy to check that commute time defined by (125]) satisfies all distance axioms. 



4-5. The Gram matrix of graph nodes 

Each vertex % G V of the graph G has an image in PM^ -1 determined by the vector 

(26) 



Vi 1 Vi 1 ~ Vn) 

The Gram matrix defined on the set of all vectors {vj}, % G V, is given by 



(27) 

The diagonal elements Gu are the first-passage times [H], the expected numbers of 
steps required for random walkers to reach the node i for the first time starting from 
any node I G V randomly chosen among all nodes of the graph G with probability 7^. 
The elements Gij, i j, estimate the expected overlap of random walks towards the 
nodes % and j starting from a node I G V randomly chosen among all nodes of the graph 
G with probability 

Discovering important nodes and quantifying differences between them in a graph is 
not easy, since the graph, in general, does not possess the structure of Euclidean space. 
The first passage time, Gu, can be directly used in order to characterize the level of 
accessibility of the node i in the graph G. Various properties of first-passage times have 
been recently studied in concern with the traffic flow forecasting [21j . in order to model 
the wireless terminal movements in a cellular wireless network [22], in a statistical test 
for the presence of a random walk component in the repeat sales price models in house 
prices [23], in the growth modelling of urban agglomerations [24] . and in many other 
works where random walks have been considered directly on city plans and physical 
landscapes. 

In contrast to all previous studies, in our approach, we use first-passage times of 
discrete time random walks in order to investigate the configuration of urban places 
represented by means of the spatial graph. 



5. Linear automorphisms of urban environments. A case of study 
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5. 1 . Spectra of cities 

If we take many, many random numbers from an interval of all real numbers 
symmetric with respect to a unit and calculate the sample mean in each case, then 
the distribution of these sample means will be approximately normal in shape and 
centered at 1 provided the size of samples was large. The probability density function 
of a normal distribution forms a symmetrical bell-shaped curve highest at the mean 
value indicating that in a random selection of the numbers around the mean (1) have 
a higher probability of being selected than those far away from the mean. Maximizing 
information entropy among all distributions with known mean and variance, the normal 
distribution arises in many areas of statistics. 




Figure 2. The probability-probability plot of the normal distribution (on the 
horizontal axis) against the empirical distribution of eigenvalues in the spectra 
of the city canal networks in Venice and Amsterdam. The diagonal line y = x 
is set for a reference. 

It is interesting to compare the empirical distributions of eigenvalues of the 
normalized Laplace operator (121 j) defined on the spatial graphs of compact urban 
patterns - the spectra of cities - with the normal distribution centered at 1. In 
Fig. [21 we have shown a probability-probability plot of the normal distribution (on 
the horizontal axis) against the empirical distribution of eigenvalues in the city spectra 
(the normal plot) of the city canal networks in Venice (96 canals) and Amsterdam (57 
canals). A random sample of the normal distribution, having size equal to the number 
of eigenvalues in the spectrum has been be generated, sorted ascendingly, and plotted 
against the response of the empirical distribution of city eigenvalues. The spectra of 
canal maintained in the compact urban patterns of Venice and Amsterdam look also 
amazingly alike and are obviously tied to the normal distribution, although these canals 
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had been founded in in the dissimilar geographical regions and for the different purposes. 
While the Venetian canals mostly serve the function of transportation routs between 
the distinct districts of the gradually growing naval capital of the Mediterranean region, 
the concentric web of Amsterdam gratchen had been built in order to defend the city. 

It is remarkable that the spectral density distributions shown in Fig. [2] are 
dramatically dissimilar to those reported for the random graphs of Erdos and Renyi 
studied by [25j I2E] . The classical Wigner semicircle distribution arises as the limiting 
distribution of eigenvalues of many random symmetric matrices as the size of the matrix 
approaches infinity, [27]. The eigenvalues of the normalized Laplace operator in a 
random scale-free graph also follow the semicircle law [28J. City spectra reveal the 
profound structural dissimilarity between urban networks and networks of other types 
studied before. 

5.2. First-passage times to ghettos 

The phenomenon of clustering of minorities, especially that of newly arrived 
immigrants, is well documented [29J (the reference appears in [30]). Clustering is 
considering to be beneficial for mutual support and for the sustenance of cultural and 
religious activities. At the same time, clustering and the subsequent physical segregation 
of minority groups would cause their economic marginalization. The spatial analysis of 
the immigrant quarters [30] and the study of London's changes over 100 years [31] 
shows that they were significantly more segregated from the neighboring areas, in 
particular, the number of street turning away from the quarters to the city centers 
were found to be less than in the other inner-city areas being usually socially barricaded 
by railways, canals and industries. It has been suggested [32] that space structure and 
its impact on movement are critical to the link between the built environment and its 
social functioning. Spatial structures creating a local situation in which there is no 
relation between movements inside the spatial pattern and outside it and the lack of 
natural space occupancy become associated with the social misuse of the structurally 
abandoned spaces. 

We have analyzed the first-passage times to individual canals in the spatial graph 
of the canal network in Venice. The distribution of numbers of canals over the range 
of the first-passage time values is represented by a histogram shown in Fig. [3jleft. The 
height of each bar in the histogram is proportional to the number of canals in the canal 
network of Venice for which the first-passage times fall into the disjoint intervals (known 
as bins). Not surprisingly, the Grand Canal, the giant Giudecca Canal and the Venetian 
lagoon are the most connected. In contrast, the Venetian Ghetto (see Fig. [3jright) - 
jumped out as by far the most isolated, despite being apparently well connected to the 
rest of the city - on average, it took 300 random steps to reach, far more than the 
average of 100 steps for other places in Venice. 

The Ghetto was created in March 1516 to separate Jews from the Christian majority 
of Venice. It persisted until 1797, when Napoleon conquered the city and demolished 
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Steps 

Figure 3. The Venetian Ghetto jumped out as by far the most isolated, despite 
being apparently well connected to the rest of the city. 



the Ghetto's gates. Now it is abandoned. 

5.3. Random walks estimate land value in Manhattan 

The notion of isolation acquires the statistical interpretation by means of random 
walks. The first-passage times in the city vary strongly from location to location. Those 
places characterized by the shortest first-passage times are easy to reach while very many 
random steps would be required in order to get into a statistically isolated site. 

Being a global characteristic of a node in the graph, the first-passage time assigns 
absolute scores to all nodes based on the probability of paths they provide for random 
walkers. The first-passage time can therefore be considered as a natural statistical 
centrality measure of the vertex within the graph, [20] . 




Figure 4. Isolation map of Manhattan. Isolation is measured by first-passage 
times to the places. Darker color corresponds to longer first-passage times. 

A visual pattern displayed on Fig. H] represents the pattern of structural isolation 
(quantified by the first-passage times) in Manhattan (darker color corresponds to longer 
first-passage times). It is interesting to note that the spatial distribution of isolation in 
the urban pattern of Manhattan (Fig. H]) shows a qualitative agreement with the map 
of the tax assessment value of the land in Manhattan reported by B. Rankin (2006) in 
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the framework of the RADICAL CARTOGRAPHY project being practically a negative 
image of that. 

Recently, we have discussed in [20] that distributions of various social variables (such 
as the mean household income and prison expenditures in different zip code areas) may 
demonstrate the striking spatial patterns which can be analyzed by means of random 
walks. In the present work, we analyze the spatial distribution of the tax assessment 
rate (TAR) in Manhattan. 

The assessment tax relies upon a special enhancement made up of the land or site 
value and differs from the market value estimating a relative wealth of the place within 
the city commonly refereed to as the 'unearned' increment of land use, [33]. The rate of 
appreciation in value of land is affected by a variety of conditions, for example it may 
depend upon other property in the same locality, will be due to a legitimate demand 
for a site, and for occupancy and height of a building upon it. 
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Figure 5. Tax assessment rate (TAR) of places in Manhattan ($ / fit 2 ) is shown 
in the logarithmic scale vs. the first-passage times (FPT) to them. 

The current tax assessment system enacted in 1981 in the city of New York 
classifies all real estate parcels into four classes subjected to the different tax rates 
set by the legislature: (i) primarily residential condominiums; (ii) other residential 
property; (iii) real estate of utility corporations and special franchise properties; (iv) 
all other properties, such as stores, warehouses, hotels, etc. However, the scarcity of 
physical space in the compact urban pattern on the island of Manhattan will naturally 
set some increase of value on all desirably located land as being a restricted commodity. 
Furthermore, regulatory constrains on housing supply exerted on housing prices by the 
state and the city in the form of 'zoning taxes' are responsible for converting the property 
tax system in a complicated mess of interlocking influences and for much of the high 
cost of housing in Manhattan, [34]. 
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Being intrigued with the likeness of the tax assessment map and the map of isolation 
in Manhattan, we have mapped the TAR figures publicly available through the Office of 
the Surveyor at the Manhattan Business Center onto the data on first-passage times to 
the corresponding places. The resulting plot is shown in Fig. [5j in the logarithmic scale. 
The data presented in Fig. |5] positively relates the geographic accessibility of places 
in Manhattan with their 'unearned increments' estimated by means of the increasing 
burden of taxation. The inverse linear pattern dominating the data is best fitted by 
the simple hyperbolic relation between the tax assessment rate (TAR) and the value of 
first-passage time (FPT), 

TAR oc (28) 
in which c ~ 120, 000 $ x Step/fit 2 is a fitting constant. 
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